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Abstract. The control of multiferroic domains through external electric fields has 
been studied by dielectric measurements and by polarized neutron diffraction on single¬ 
crystalline TbMnOa- Full hysteresis cycles were recorded by varying an external field 
of the order of several kV/mm and by recording the chiral magnetic scattering as well 
as the charge in a sample capacitor. Both methods yield comparable coercive fields 
that increase upon cooling. 
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Control of multiferroic domains by external electric fields in TbMnOs 2 

1. Introduction: chiral domains in multiferroic compounds 

Multiferroics continue attracting strong interest due to their application potential in data 
storage and processing DEI. In particular the more recently discovered multiferroic 
materials, in which the ferroelectric polarization is directly coupled to the magnetic 
order |2], are promising, because they allow one to control the magnetic structure 
through an external electric held. In these materials the magnetic structure is rather 
complex with non-collinear and frequently chiral arrangement of magnetic moments. 
Such a non-collinear arrangement can induce a local electric dipole through the inverse 
Dzyaloshinski-Moriya mechanism, P oc fij x [Si x Sj), where Cij is the connecting vector 
between two neighboring spins Si and Sj [3]. For example in a cycloid these local 
dipoles add up to macroscopic ferroelectric polarization. This mechanism perfectly 
describes the multiferroic phases in many of the transition-metal oxides multiferroics 
such as TbMnOsIl], Ni3V208 |5] and MnW 04 |6]. The inverse Dzyaloshinski-Moriya 
interaction is not the only possibility to induce ferroelectric polarization through a 
canted spin arrangement mi- For example for lower symmetry such as in monoclinic 
systems polarization may also appear parallel to S'* x Sj, which applies to the multiferroic 
phases in CuFe 02 [Tl E], CuCr02[in], RbFe(Mo 04 ) 2 [II] and NaFeSi 206 [T 2 ] . 

In the multiferroic phase there are at least two different electric domains 
corresponding to the direction of the ferroelectric polarization. These electric domains 
are coupled with the magnetic structure, which is chiral in most cases. A cycloid does 
not exhibit a scalar chirality but a vector chirality, which can be directly associated with 
the ferroelectric polarization. Therefore, the magnetoelectric coupling in the multiferroic 
phase allows one to control the chiral magnetism with the aid of an external electric field. 
In a future application one can add an exchange bias coupled ferromagnetic layer on 
top of the multiferroic compound in order to read the multiferroic state [13]. Studying 
the multiferroic domains in a bulk material is, however, more difficult and requires a 
microscopic technique such as second harmonic generation miis] or polarized neutron 
diffraction [TB] . 

Polarized neutron diffraction is ideally suited to study chiral magnetism [T6|. This 
can be easily understood, because the neutron is a chiral object itself when its magnetic 
moment and momentum are parallel. In general magnetic neutron diffraction measures 
the Fourier component of the magnetization distribution, M, but only the components 
perpendicular to the scattering vector, Q, contribute, Mj_ with |M_lP = \My\‘^ + 

Here we use the common coordinate system in polarized neutron scattering experiments 
—X II Q, T vertical to the scattering plane, and y = zx x. The chiral magnetic scattering 
is defined as = —i[M± x as it can be illustrated by considering an ideal 

helix. For such ideal helix and Q parallel to the propagation vector is maximal and 
of the same size as |Mxp. Using spherical polarization analysis one may analyze the 
full polarization matrix, i.e. analyzing the neutron polarization in the outgoing beam 
parallel to x, y, or z for any polarization in the incoming beam. The intensity in such 
a channel, auv with u, v = x, y, or z, can be recorded for up and down polarization. 
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respectively, resulting in four values in each channel. We denote the down polarization 
by an overbar. For the longitudinal cases, CFyy and one may distinguish spin- 
flip (SF) and non-spin-flip (NSF) processes. The chiral components can be directly 
determined in various channels, most easily they are obtained in the xx channels where 
all the magnetic scattering contributes to the SF scattering. For these two SF intensities 
one hnds Oxx = M± ■ — i{M± x Mf)x and axx = M± ■ Mf_ + i{M± x Mf)x- In one 

case the chiral contribution adds to while it is subtracted in the other case. In 

the case of the ideal helix with the propagation vector parallel to Q the total intensity 
will thus appear only in one of the two SF xx channels. However, this applies only to 
a mono-domain sample, while for equal domain distribution the two axx and a^x SF 
channels become again identical. Comparing these two SF channels gives thus a direct 
method to determine the domain distribution in a chiral magnet. [ini [HI EHl DS]. 

Concerning the multiferroic materials neutron polarization analysis was hrst applied 
for TbMnOs to show that one can control the chiral domains by cooling in an external 
electric held ng. Similar temperature dependent measurements were performed for 
LiCu202 [20] and MnW 04 [HI |2T]. For MnW 04 neutron polarization analysis was 
also used to follow multiferroic hysteresis cycles by varying the external electric held at 
constant temperature and recording the chiral domain distribution [T8l l2^ 123] . 

TbMnOs is the prototypical material for the new class of chiral multiferroic 
materials, because it exhibits sizable ferroelectric polarization and large magnetoelectric 
coupling combined with a moderately complex crystal and magnetic structure [U |2Tl |25] . 
Magnetic order in TbMnOs sets in at Tn=42K leading to a longitudinal spin-density 
wave (SDW) with a propagation vector of kinc={fi ~0.28 0) in reduced lattice units of the 
Pbnm structure. Upon further cooling below Tc=27.3K a second magnetic transition 
occurs into an elliptical cycloid phase with nearly the same propagation vector. In the 
SDW phase spins are aligned parallel to the b direction, while they rotate in the be plane 
in the cycloidal phase. This second transition is accompanied by the onset of ferroelectric 
order, which is perfectly explained by the inverse Dzyaloshinski-Moriya mechanism. The 
same mechanism also describes the rotation of the ferroelectric polarization in external 
magnetic helds [26] , as the cycloid plane flops simultaneously from the be to the ab plane 
[27]. The multiferroic hysteresis loops have not been examined in detail for TbMnOs, 
only a single loop was studied by second harmonic generation on a thin him |28j . 

Here we study the control of the multiferroic domain population by application of 
an external electric held in TbMnOs. We directly measure the ferroelectric polarization 
by counting the surface charge in a capacitor and we determine the chiral domain 
distribution by polarized neutron dihraction. Both methods document that it is possible 
to switch the multiferroic domains in bulk TbMnOs by varying the external held at 
constant temperature, but only sufficiently close to the ferroelectric transition. 
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Waveformgenerator High Voltage Amplifier Sample and Electrometer 

Nl 6251 BNC Trek 609E reference capacitor Keithley 6517b 


A/D Converter 
Nl 6251 BNC 


Figure 1. Scheme of the setup used for dielectric measurements of the multiferroic 
domains in TbMnOa. It consists of a modified Sawyer-Tower circuit. 


2. Experimental 

A large single crystal of TbMnOs was grown with the travelling solvent floating-zone 
technique in a mirror furnace (Crystal Systems Incorporated FZ-T-IOOOO-H-VI-VP). 
The magnetic transitions were determined by SQUID measurements on two identically 
grown crystals to occur at TAr=42.0(2) and T(7=27.3(2) K in good agreement with 
literature I11121]. All crystals used in the neutron diffraction and macroscopic studies 
in this work were cut from the same single crystal. 

Polarized neutron diffraction experiments were performed with the polarized single¬ 
crystal diffractometer POLI-HEIDI installed at the hot source of the Forschungsreaktor 
Miinchen II (FRM-II). The neutron polarization and analysis is accomplished by ^He- 
spin hlters in the incoming and scattered beams [221 EDI El]- In order to control the 
decay of the hrst hlter polarization the incoming beam polarization was recorded by a 
transmission monitor. The polarization analysis of the scattered beam was verihed by 
analyzing a purely nuclear Bragg reflection. The direction of the neutron polarization 
at the sample was controlled with the CRYOPAD setup [16] . 

In order to study the effect of crystal thickness and to apply an electrical held 
parallel to the c direction, two plate-shaped single-crystals were cut perpendicular to the 
c axis with a thickness of 3.1 mm and 0.9 mm, respectively. The electric held was applied 
with a capacitor build from aluminum plates and nylon screws. The sample was mounted 
in the (2 0 l)/(0 1 0) scattering plane. The strong magnetic rehection Q = [2 0.28 1) was 
used for determining the chiral ratio. Note that a large h component of the magnetic 
Bragg rehection is needed in this experiment, because both components of the be cycloid 
must strongly contribute to the scattering, which requires that the b and c axes are nearly 
perpendicular to Q. The intensity of the thin sample at this rehection was rather low, 
50cts/s. The sample was cooled with a standard FRM-II type top-loading closed-cycle 
refrigerator, which stabilizes temperatures within 0.1 K. The challenge of measurements 
with high electric helds in croygenic environments consists in hne-tuning the pressure 
of helium exchange gas continuously on a thin line where it is possible to apply several 
kV/mm without electric discharge while the sample temperature is still controllable. 
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Figure 2. (Color online) The magnetic scattering at the (2 0.28 1) Bragg reflection 
in TbMnOa was recorded in the two SF channels for neutron polarization along 
the X direction. During cooling through the magnetic transitions, electric fields of 
±i.3kV/mm were applied in (a) and (c), while the data in (b) were taken in zero field. 
The intensities of the two SF channels axx and Uxx interchange upon reversing the 
field, which indicates the population of opposite chiral domains. 
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Much efforts were needed to determine the breakdown pressure versus voltage curves 
before starting the experiment. A temperature offset of about 1 K was observed for the 
phase transitions measured with applied field and those reported in the literature. This 
offset arises from the insufficient thermal contact between the temperature sensor and 
the isolating sample, when the exchange gas is lowered to allow for the application of 
electric field. In the following we assume this temperature offset to be constant in the 
temperature range from 10 to 50 K and correct the data correspondingly. The chiral 
ratio was determined by measuring the two SF axx and a^x channels as described above. 

The macroscopic ferroelectric polarization was measured with a modified Sawyer- 
Tower circuit, in which the surface charge is determined indirectly with an electrometer 
(Keithley 6517b), see Figure 1. A thin plate of 0.2 mm thickness and a contacted 
surface area of 9 mm^ was cut from the same large crystal as those used in the neutron 
diffraction experiments. The crystal was contacted on both sides with conductive silver 
paint in order to form the sample capacitor, which is studied in respect to the reference 
capacitor in the Sawyer-Tower circuit. We applied electric fields of up to 3 kV/mm 
and recorded the hysteresis loops in the quasistatic configuration. Running a full loop 
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Figure 3. (Color online) Chiral ratio measured at T=25K after cooling the thick 
sample in small electric fields from the paramagnetic (50 K) into the chiral magnetic 
phase (25 K). The sequence was measured with reducing and alternating field strengths. 
At low electric fields the chiral ratio determined at 25 K shows a linear dependence on 
the electric field applied during cooling. 

corresponds to one second. The high voltage was generated by symmetric triangular 
waveforms which have been amplihed by a factor of 1000 in a Trek 609E high voltage 
ampliher. In a vacuum of about 10“® mbar it was possible to apply electric helds of the 
order of 5kV/mm. The results are limited by hnite conductivity losses, which turned 
out to be small. The sample capacitor was mounted on the cold huger of a closed cycle 
refrigerator, which allows to reach base temperatures of 5 K. The control of temperature 
and the runs to record the hysteresis loops were fully automatized. Comparing several 
hysteresis loops no ongoing fatigue effects could be detected, but fatigue effects for an 
asgrown crystal cannot be excluded. Therefore, it is possible to average 4 hysteresis 
loops at each temperature. This procedure effectively suppresses noise which partially 
stems from the triboelectric effect. 

3. Results 

3.1. Polarized neutron diffraction experiments 

In a preliminarily test the sample was cooled from 50 K in positive and negative electric 
helds of ±1.3kV/mm and in zero held, see Figure 2. The sample develops a clear 
preference for one chiral domain depending on the sign of the electric held. The intensity 
of the two spin-hip channels is interchanged by reversing the electric held. The chiral 
ratio, Tchir = , was calculated from the data for cooling in positive held and 

roughly follows the development of the electric polarization as it has been previously 
reported by the same type of experiment [12]. Note that even for a mono-domain 
state Tchir will deviate from ±1 due to the hnite geometry factors in TbMnOa, which 
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Figure 4. (Color online) Hysteresis loops in TbMnOa obtained by measuring the 
chiral ratio as a function of external electric field at different temperatures for two 
different samples at POLI-HEIDI (thicknesses of 0.9 and 3.1mm). Open triangles 
denote the data taken immediately after cooling and filled symbols denote the data 
taken to characterize the full hysteresis loops. Triangles plotted in orange show the 
three loops taken with the thicker crystal. Crosses denote the coercive field values. 
The loops were recorded after zero-field cooling from 50 K. 


furthermore is not a perfect cycloid [2S]. With the elliptical magnetic structure we 
calculate an ideal mono-domain chiral ratio rchir=±0.93 for Q = [2 0.28 1). 

The effect of small electric helds applied during cooling from the paramagnetic into 
the chiral magnetic phase is depicted in Figure 3. The sequence was measured with 
reducing and alternating held strengths. Very small helds of a few tenths of V/mm 
applied upon cooling are sufficient to produce hnite preference for one chiral domain. 
However, saturation is not reached for such small helds. Instead, the chiral ratio shows 
a linear dependence on the electric held applied during the cooling, at least in this low- 
held range. Note, however, that the chiral response of the samples is asymmetric with 
respect to the applied held, as it even exhibits hnite values for zero electric held. 

In the following we focus on the possibility to reverse chiral magnetic domains 
by varying the electric held at constant temperature. This represents the multiferroic 
hysteresis loops, i.e. the control of magnetic order by an external electric held. The 
obtained hysteresis loops are shown in Figure 4. All hysteresis loops were recorded after 
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Figure 5. (Color online) (a) Temperature dependence of the remanent ferroelectric 
polarization in TbMnOa as obtained from the dielectric hysteresis loops shown in 
Fig. 6. In comparison we also show the temperature dependence of the remanent 
chiral ratios observed in polarized neutron diffraction. In both parts larger triangles 
correspond to the neutron results, and the tip of the triangle indicates the direction of 
the applied electric fields, (b) Temperature dependence of the coercive fields obtained 
in the neutron diffraction and dielectric measurements. The positive and negative 
coercive fields obtained from the chiral-ratio loops match perfectly underlining the 
symmetric character of the sample aside from the results obtained for the thicker crystal 
with neutron diffraction. Furthermore, the coercive fields obtained in the neutron 
experiment do not depend on sample thickness. The coercive field obtained from the 
dielectric experiments are considerably higher. 


zero-field cooling from 50 K to the respective temperature. To ensure that there was no 
electric held applied during cooling the sample surfaces were short-circuited for most of 
the loops shown in Figure 4. Only the curve taken at 19 K was obtained after cooling the 
sample with open contacts, which seems to induce a different direction of preference. It 
is a remarkable hnding that the initial curves always start at hnite values even for cooling 
with short-circuited surface contacts. The hysteresis loops were recorded by varying the 
electric held in hnite steps and by recording the neutron intensities at constant helds. 
Due to the limited count rates, it was necessary to count several minutes per point. 
Therefore, the shown hysteresis loops are quasi-static for relaxation times of the order 
of minutes and below. 

Most recorded hysteresis loops show a symmetric shape. Both the coercive helds 
and the saturation values are in good accordance between the two directions. The 
coercive held increases with decreasing temperature while an approximately constant 
saturation value of rchir can be reached at low temperatures by applying a high enough 
held, see Fig. 5(a) and (b). Note, however, that the chiral ratio describes the domain 
population multiplied with the intrinsic chiral fraction of the magnetic scattering and, 
therefore, it is not expected to continuously increase upon cooling while the ferroelectric 
polarization follows Si x Sj. The domain dynamics of the sample becomes stiher when 
cooling deeper into the multiferroic phase. This behavior qualitatively resembles the 
results reported for MnW 04 imiiH]. At 9 K the domains are strongly pinned; therefore. 
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the electric field causes an electrical breakdown before the chiral ratio can be affected. 
Already at 14 K the maximal applicable electric field hardly surpasses the coercive 
field needed to start reverting the multiferroic domains within the time scales of the 
experiment. 

The hysteresis loops of the two samples of different thicknesses match well each 
other concerning the coercive fields aside from a more pronounced asymmetry for the 
thicker crystal. In order to apply the same electric field strength to the thicker sample a 
larger voltage has to be applied. For that reason it is not possible to reach as high electric 
fields for the thicker sample as for the thinner sample. Therefore, we registered only 
three hysteresis loops at higher temperatures for the sample, which was 3.1mm thick. 
The coercive fields for all temperatures are shown in Figure 5 (b). The positive and 
negative coercive fields match well underlining the symmetric character of the thinner 
samples. The coercive field shows roughly a linear dependence on the temperature. 
There is less agreement concerning the remanent values of the chiral ratio observed in 
the two samples of different thickness. The thicker sample of the neutron experiment 
remains less polarized. 

3.2. Ferroelectric hysteresis loops detected by capacitor measurements 

The ferroelectric polarization of TbMnOs was studied macroscopically with a modified 
Sawyer-Tower circuit. The ferroelectric polarization P versus electric field E hysteresis 
loops obtained by averaging four successive cycles are shown in Fig. 6 for different 
temperatures. The hysteresis loops are not fully fiat in their saturation region due to 
the dielectric background contribution, as it is usually observed in ferroelectrics, and 
due to finite conductivity effects. Therefore we have substracted a loop taken in the 
paraelectric phase from all the other loops. Finite slopes then still arise from the minor 
differences in conductivity and dielectric constant occurring with the variation of the 
temperature. For clarity we show only the data taken in 0.2 K temperature steps, but 
intermediate temperatures were also measured. 

The hysteresis loops in Fig. 6 show that the ferroelectric transition sets in at 
27.3(1) K in perfect agreement with the susceptibility data and literature. Below this 
temperature we find remanent polarization as well as finite width and area in the loops. 
The slight change in the slope of the entire loop can be attributed to the maximum of 
the dielectric constant at Tc |1]. 

The area of the loop considerably increases upon cooling. The remanent 
polarization is almost identical to the saturation polarization after correcting for the 
finite slope. The temperature dependence of the remanent polarization matches the 
temperature dependence of the polarization measured by the pyrocurrent method |H |21] 
as long as temperatures stay above 24 K, see Fig. 5 and 6. This indicates that for 
temperatures larger than 24 K the polarization can be fully controlled with the maximum 
applied electric field of 3 kV/mm and within a time of about one second. In contrast the 
data taken below 24 K show that saturation polarization cannot be reached indicating 
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Figure 6. (Color online) Dielectric hysteresis loops recorded by a modified Sawyer- 
Tower circuit on a 0.2 mm thick sample. The colors indicate the temperature 
varied within each panel in 0.2 K steps from the colder (blue) to the warmer (red) 
temperatures. A linear background determined above the ferroelectric transition was 
substracted. 

that the relaxation times surpass the times of the measurement. Between 23 and 24 K 
one may recognize that saturation is not reached before the held starts to decrease from 
its maximal values, and below 23 K the polarization continues to increase in absolute 
size even when the electric held is already decreasing. This clearly shows that the 
measurement was too fast to allow the TbMnOa domains to fully relax. Note that the 
neutron experiment due to the low count rates was sufficiently slow to allow for full 
relaxation to lower temperatures. 

The temperature dependence of the coercive held is shown in Fig. 5(b). In this 
dielectric experiment there is no signihcant diherence in the coercive helds obtained 
for the two held directions, which agrees to the neutron experiments and also to the 
symmetric values of the remanent polarization. In the studied temperature interval the 
coercive held increases linearly with temperature after a small jump at the transition. 
Compared to the neutron study the coercive helds are considerably larger in spite of 
the fact that all the crystals were cut from the same large single crystal. However, 
the much thinner sample used in the dielectric study can more strongly suher from 
the cutting process. In addition fatigue ehects can be stronger in the sample used for 
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the macroscopic measurement as it went through a much larger number of thermal 
and electrical cycles before the shown experiments were performed. Coercive helds in 
TbMnOs are seem to be sample dependent. 

4. Conclusions 

In conclusion we have shown by polarized neutron diffraction and by macroscopic 
measurements that the population of multiferroic domains in bulk TbMnOs can be 
controlled through a moderate electric held of the order of several kV/mm. Neutron 
diffraction determines the reversion of the chiral magnetism while the macroscopic 
measurement senses the ferroelectric polarization. The coercive helds in both studies 
are of the same order but signihcant diherences can be ascribed to diherences in the 
samples used. The chiral magnetic and polarization hysteresis loops conhrm once more 
the close coupling of the two properties. 
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